A re-annotation pipeline for Illumina BeadArrays: improving the interpretation of gene expression data.

Illumina BeadArrays are among the most popular and reliable platforms for gene expression profiling. However, little external scrutiny has been given to the design, selection and annotation of BeadArray probes, which is a fundamental issue in data quality and interpretation. Here we present a pipeline for the complete genomic and transcriptomic re-annotation of Illumina probe sequences, also applicable to other platforms, with its output available through a Web interface and incorporated into Bioconductor packages. We have identified several problems with the design of individual probes and we show the benefits of probe re-annotation on the analysis of BeadArray gene expression data sets. We discuss the importance of aspects such as probe coverage of individual transcripts, alternative messenger RNA splicing, single-nucleotide polymorphisms, repeat sequences, RNA degradation biases and probes targeting genomic regions with no known transcription. We conclude that many of the Illumina probes have unreliable original annotation and that our re-annotation allows analyses to focus on the good quality probes, which form the majority, and also to expand the scope of biological information that can be extracted.

Autophagy mediates the mitotic senescence transition.

As a stress response, senescence is a dynamic process involving multiple effector mechanisms whose combination determines the phenotypic quality. Here we identify autophagy as a new effector mechanism of senescence. Autophagy is activated during senescence and its activation is correlated with negative feedback in the PI3K-mammalian target of rapamycin (mTOR) pathway. A subset of autophagy-related genes are up-regulated during senescence: Overexpression of one of those genes, ULK3, induces autophagy and senescence. Furthermore, inhibition of autophagy delays the senescence phenotype, including senescence-associated secretion. Our data suggest that autophagy, and its consequent protein turnover, mediate the acquisition of the senescence phenotype.

Detecting the coevolution of biosequences--an example of RNA interaction prediction.

A probabilistic graphical model is proposed in order to detect the coevolution between different sites in biological sequences. The model extends the continuous-time Markov process of sequence substitution for single nucleic or amino acids and imposes general constraints regarding simultaneous changes on the substitution rate matrix. Given a multiple sequence alignment for each molecule of interest and a phylogenetic tree, the model can predict potential interactions within or between nucleic acids and proteins. Initial validation of the model is carried out using tRNA and 16S rRNA sequence data. The model accurately identifies the secondary interactions of tRNA as well as several known tertiary interactions. In addition, results on 16S rRNA data indicate this general and simple coevolutionary model outperforms several other parametric and nonparametric methods in predicting secondary interactions. Furthermore, the majority of the putative predictions exhibit either direct contact or proximity of the nucleotide pairs in the 3-dimensional structure of the Thermus thermophilus ribosomal small subunit. The results on RNA data suggest a general model of coevolution might be applied to other types of interactions between protein, DNA, and RNA molecules.